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Abstract 

We present a model of financial markets originally proposed for a turbulent flow, as 
a dynamic basis of its intermittent behavior. Time evolution of the price change is 
assumed to be described by Brownian motion in a power-law potential, where the 
'temperature' fluctuates slowly. The model generally yields a fat-tailed distribution 
of the price change. Specifically a Tsallis distribution is obtained if the inverse tem- 
perature is ^-distributed, which qualitatively agrees with intraday data of foreign 
exchange market. The so-called 'volatility', a quantity indicating the risk or activity 
in financial markets, corresponds to the temperature of markets and its fluctuation 
leads to intermittency. 
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1 Introduction 



Financial returns are known to be non-gaussian and exhibit fat-tailed distri- 
bution ([l|)-(|T5D. The fat tail relates to intermittency — an unexpected high 
probability of large price changes, which is of utmost importance for risk anal- 
ysis. The recent development of high-frequency data bases makes it possible 
to study the intermittent market dynamics on time scales of less than a day 
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Using foreign exchange (FX) intraday data, Miiller et al. (|P2|) showed that 
there is a net flow of information from long to short timescales, i.e., the be- 
havior of long-term traders influences the behavior of short-term traders. Moti- 
vated by this hierarchical structure, Ghashghaie et al. (f|) have discussed analo- 



gies between the market dynamics and hydrodynamic turbulence (IB; |T7|), and 
claimed that the information cascade in time hierarchy exists in a FX mar- 
ket, which corresponds to the energy cascade in space hierarchy in a three- 
dimensional turbulent flow. These studies have stimulated further investiga- 
tions on similarities and differences in statistical properties of the fluctuations 
in the economic data and turbulence (^)-(|H]). Differences have also emerged. 
Mantegra and Stanley (|| 0) and Arneodo et al. (|) pointed out that the time 
evolution, or equivalently the power spectrum is different for the price differ- 
ence (nearly white spectrum) and the velocity difference (J 1 / 3 spectrum, i.e., 
— 5/3-law for the spectrum of the velocity). Moreover, from a parallel analysis 
of the price change data with time delay and the velocity difference data with 
time delay (equivalent to the velocity difference data with spatial separation 
under the Taylor hypothesis (|16D ), it was shown that the time evolution of the 
second moment and the shape of the probability density function (PDF), i.e., 
the deviation from gaussian PDF are different in these two stochastic processes 
((?]). On the other hand, non-gaussian character in fully developed turbulence 
(|T6|) has been linked with the nonextensive statistical physics (IT 



As dynamical foundation of nonextensive statistics, Beck recently proposed 
a new model describing hydrodynamic turbulence ( PT[ |2"2"| ). The velocity dif- 
ference Av of two points in a turbulent flow with the spatial separation Ar 
is described by Brownian motion (an overdamped Langevin equation (H)) 
in a power-law potential. Assuming a ^-distribution for the inverse temper- 
ature, he obtained a Tsallis distribution fll8|) for Av. However, if we take 
into account the almost uncorrelated behavior of the price change ©-(0), 
the picture by means of the Brownian motion seems to be more appropriate 
for the market data rather than turbulence. Moreover, the description by the 
Langevin equation is able to relate the PDF of the price change to that of the 
volatility, a quantity known as a measure of the risk in the market. Thus we 
applied the model to FX market dynamics by employing the correspondence 
by Ghashghaie et al. (^). 



2 Model 



We substitute the FX price difference Z(t) = y(t + At) — y(t) with the time 
delay At for the velocity difference Av with the spatial separation Ar. Beck's 
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model for turbulence then reads 



f =7F(Z) +*(«), 



where 7 > is a constant, and R(t) is gaussian white noise corresponding 
to the temperature kT, satisfying (R(t)R(t')) = 27 kT5(t — f). The 'force' 
F = —dU(Z)/dZ is assumed to be obtained by a power-law potential U(Z) = 
C\Z\ 2a with an exponent la, where < a < 1 and C is a positive constant. 
That is, the system is subject to a restoring force proportional to the power 
of the price difference, \Z\ 2a besides the random force R(t). Especially when 
a = 1, the restoring force is linear to Z. Under a constant temperature kT, 
Eq. (1) leads to a stationary (i.e., thermal equilibrium) distribution of Z as 

D -pu(Z) 

^ - jhmrz - ^m^e-^, ( 2) 



where (3 = 1/kT is the inverse temperature (|2T|) . The 'local' variance of Z, 
which is defined for a fixed value of (3, is obtained from the conditional prob- 
ability in Eq. (2) as 

Z 2 )^ = J Z 2 P At (Z\(3)dZ = ^^(Cpy 1 ^ oc (kT) l/a . (3) 
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We define the volatility by the square root of the local variance of Z (see 
Eq. (9)). When a = 1 and C = 1/2, (Z 2 )^ coincides with the temperature 
and the conditional PDF reduces to gaussian, while for a / 1, {(Z 2 )^}" is 
proportional to kT. 



Let us assume that the 'temperature' of the FX market is not constant and 
fluctuates in larger time scales, and j3 is, just as in Beck's model for turbulence, 
X 2 -distributed with degree n: 

^mi^r^M-S)' " >2 - (4) 



where T is the Gamma function, /3 is the average of the fluctuating (3 and n 
relates to the relative variance of f3: 

(P)=Po, 0)2 =-. (5) 
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Equation (4) implies that the local variance (Z 2 ) ~ = v fluctuates with the 

distribution p(v) oc v~( a+1 ' fAt{v~ a ) ■ The conditional probability in Eq. (2) 
together with Eq. (4) yields a Tsallis-type distribution ( |i~8| for the ultimate 
PDF of Z: 



l i+ 2a-(q-l) J 

1 < Cf3 2a( q -l) \ 1/2a r(^) 

^ 1 2« - (g - 1) / r(i)r(^-i)' 10 

where Tsallis' nonextensivity parameter g is defined by 

9 = 1 + -^-, (8) 
an + 1 



which satisfies 1 < q < 5/3 because of < a < 1 and n > 2. Since g > 1, the 
distribution of Z exhibits power-law tails for large Z: P& t {Z) ~ ^^ 2q /(?^ 1 ) = 
2j-{an+x) _ fj ence; TTjth moment (Z m ) = J Z m P 1 \ t (Z)dZ converges only for 
m < an. In the limit of q — > 1, P& t {Z) in Eq. (6) reduces to the canonical 
distribution of extensive statistical mechanics: lim^x Pa<(Z) = Pm{Z\(3q). 



3 Results and discussions 



We have applied the present model to the same FX market data set as used 
in Ref. (f|) (provided by Olsen and Associates (^5|) which consists of 1 472 241 
bid-ask quotes for US dollar-German mark exchange rates during the period 
October 92 - September 93). The volatility is often estimated by the standard 
deviation of the price change in an appropriate time window (0). Employing 
this definition, we express the volatility in terms of the local standard deviation 
of the price change Z(ti) = y(ti + At) — y(ti) as 



ttE{^)} 2 - (9) 

iV 3=1 




Here the window size has been chosen as iV = 35. Since V~ 2a corresponds 
to (3 (see Eq. (3)) which is x 2 -distributed, n can be explicitly obtained from 
the relative variance of (3 using Eq. (5). Thus, there is only one adjustable 
parameter among (a, q, n), because we have another relation, Eq. (8). In other 
words, the PDF, P^ t (Z) of the price change and the PDF, fAt(P) of the inverse 
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power V 2a = (3 of the volatility are determined simultaneously once the value 
of a has been specified. 
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Fig. 1. Data points: standardized PDF PAt(Z) °f price changes 
Z(t) = y(t + At) - y(t) for time delays At = 320 s, 1 280 s, 5 120 s, 20 
480 s (from top to bottom). The middle prices y(t) = (yud(t) + Vask(t))/2 have been 
used (data provided by Olsen and Associates (|25|)), The PDF has been obtained in 
a similar way as Ref. (H), i.e., we have selected the complete set of non-overlapping 
records separated by a time interval At(l ± e) with the tolerance < e < 0.035. 
The number of the available data points thus decreases with At, which is 54 087, 
16 338, 4 359, 1 114 (from top to bottom). For better visibly, the curves have been 
vertically shifted with respect to each other. Full lines: theoretical expression given 
in Eq. (6) in text, where (a,q,n) is (0.918, 1.393, 4), (0.832, 1.322, 5), (0.565, 
1.101, 18), (0.540, 1.054, 35) (from top to bottom), n is explicitly decided from the 
variance of the volatility data (the closest integer satisfying Eq. (5)), q is from Eq. 
(8), so that we only need to adjust a by least square fitting. 

The PDF's with time delay At varying from five minutes up to approximately 
six hours are displayed in Fig.l together with theoretical curves obtained from 
Eq. (6). As the time scale At increases, n increases, while a and q decrease. The 
nonextensivity parameter q tends to the extensive limit: q — > 1 as At increases. 
Using the same parameter values (a,q,n), the PDF's of (3 are compared in 
Fig. 2. The average 'temperature' in the market 1/(3q ~ (Z 2 ) a ~ (At) QI > 2 
increases with At since a( 2 is positive. (We obtained the scaling exponent 
djj; ^) £ « 0.85, which is larger than 2/3 obtained for turbulence.) In contrast, 
the fluctuation of the temperature increases with decreasing At because the 
variance of the inverse temperature is proportional to 2/n. The smaller values 
of n imply the stronger intermittency which occurs in small time scales. The 
intermittent character in the price change can be seen as a fat tail of Pm(Z) 
in Fig. 1. Also in Fig. 2, the peak of fAt(P) shifts to smaller (3 as At decreases 
(from d to a in Fig. 2), which means relatively high temperatures are realized 
more frequently in short time scales. 
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Fig. 2. Rescaled probability density fAt(P) of the quantity V~ 2a = f3, where V is 
the volatility defined by the local standard deviation of the price change Z(t). (See 
the legend of Fig. 1 for the definition of Z(t).) Thin lines: obtained from the volatility 

data, which was defined by V(ti) = ^Y^f=i {^i^N(i-i)+j)} 2 ~/N with = 35 so 
that the windows for different volatility data points do not overlap. The number 
of the available data points is 1 545, 466, 124, 24 (from a to d). The data set has 
been rescaled so that the average of /3 is unity: (/?) = (V~ 2a ^ = 1. Thick lines: 
X 2 -distribution for j3 given in Eq. (4) with (3q = 1. The same parameter values as in 
Fig. 1 are used, where n is the integer closest to 2/ ((/3 2 ) — l) (see Eq. (5) in text), 
a has been adjusted by least square fitting of P^tiZ) in Fig. 1. (a,q,n) is (0.918, 
1.393, 4), (0.832, 1.322, 5), (0.565, 1.101, 18), (0.540, 1.054, 35) (from a to d). 

It should be noted that the PDF in Eq. (6) with a = 1 reduces to Student's 
t-distribution (§), which has often been used to characterize the fat tails fll5l). 
When a — 1, there is no adjustable parameter because n is decided from Eq. 
(5). The FX market system is then subject to a restoring force linear to the 
price change and the volatility is proportional to the temperature. We have 
found that the PDFs of Z and (3 reproduce, although very roughly, the Olsen 
and Associates' data points even if a is fixed at a = 1. However, adjusting 
the parameter a improved the line shape of Pm(Z) in a range close to Z = 0. 
The data points in Fig. 1 exhibit a cusp at Z = for large time scales At, 
which implies a singularity of the second derivative of the PDF at Z = 0. The 
larger reduction of a from unity leads to the stronger singularity. (Note that 
the factor \Z\ 2a ~ 2 arises from d 2 PAt{Z) / dZ 2 .) Thus the better fitting for large 
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At was obtained from a reduced value of a. The trend of decrease in a with 
increasing At was observed for the turbulent flow ( pT|) as well. However, the 
deviation from a = 1 is much smaller than the present case and the cusp is 
invisible. 

Ghashghaie et al. @ have used a model for turbulence by Castaing et al. (|17D, 
in which a log-normal distribution has been assumed for the local standard 
deviation of the price change. The present model reduces to the model by 
Ghashghaie et al. if the stochastic process given in Eq. (1) is assumed with a = 
1 (the local variance of Z is then proportional to kT) and the ^-distribution 
for (3 (the inverse of the local variance of Z) is replaced by the log-normal 
distribution for \JkT (the local standard deviation of Z). Although no analytic 
expression like Eq. (6) for P At (Z) is obtained, a similar qualitative explanation 
can be applied to their model: The volatility (or equivalently, the square root 
of the temperature) fluctuates slowly with a log-normal distribution, and the 
smaller time scale corresponds to the larger variance of the logarithm of the 
volatility. (The variance is denoted by A 2 in Ref. 4.) However, the power-law 
behavior of the tail of the volatility distribution (|TJ]) can be better described 
by the x 2 -distribution for the inverse of the variance. 

We have proposed the stochastic process described by Eq. (1) for FX mar- 
ket dynamics in small time scales. In fact, Eq. (1) is the simplest stochastic 
process which can realize the thermal equilibrium distribution, Eq. (2) in the 
power-law potential. More realistic and more complicated processes that as- 
sure convergence to Eq. (2) at local temperatures might be possible. Mantegna 
and Stanley have proposed a different stochastic model of the price change, 
which is described by a truncated Levy flight (TLF) flT3| ; ^). The model, yield- 
ing approximately a stable distribution, well reproduces the self-similar prop- 
erty of PDF at different time scales: At = 1 min to 1 000 min. However, 
the parameters a and 7 characterizing the stable distribution fluctuate for 
larger (monthly) time scales @, where 7 gives a measure of the volatility. 
In other words, the ultimate distribution (let us denote it Pj^ ¥ (Z)) should 
be obtained, like Eq. (6), from the weighted average over these parameters. 
A difference between P]^ F (Z) and P/\ t {Z) is that P]^ F (Z) has no cusp at 
Z = 0: Plf F (Z) ~ 1 - AZ 2 for \Z\ < 1, whereas P At (Z) ~ 1 - A\Z\ 2a } and 
the present FX date set indeed exhibits a cusp as seen in Fig. 1. 

Finally, a fundamental question has been left open: How to derive theoretically 
the increasing trend of fat-tailed character with decreasing time scale, i.e., 
At-dependence of the parameters (a, q, n). An attempt to derive the non- 
gaussian fat-tailed character in small time scales was made by Friedlich et al. 
recently @. They derived a multiplicative Langevin equation from a Fokker- 
Planck equation and showed that the equation becomes more multiplicative 
and hence fat-tailed as At decreases. Clarifying the relation between their 
multiplicative Langevin equation and Eq. (1) (the latter being rather simple 
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although including the fluctuating temperature) should be important as well 
as interesting. 
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